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RANDOMNESS and COLLECTIVITY in NUCLEAR 
STRUCTURE: 

Three theoretical puzzles 


Vladimir ZELEVINSKY and Alexander VOLYA 
Department of Physics and Astronomy and 
National Superconducting Cyclotron Laboratory, 

East Lansing, MI f8824-1321 USA 

We show and interpret three examples of nontrivial results obtained in numerical 
simulations of many-body systems: exponential convergence of low-lying energy 
eigenvalues in the process of progressive truncation of huge shell-model matrices, 
apparently ordered spectra generated by random interactions, and regular behav¬ 
ior of complex many-body energies in a system with single-particle orbitals in 
continuum. The possible practical applications and new approaches are suggested. 


1 Introduction 

As we discussed§ at the previous conference (S. Agata sui due Golfi, 1998), the 
ideas of quantum chaos significantly advance our understanding of many-body 
quantum systems. The new stage of this development is related to the role of 
incoherent interactions between the constituents. Apart from great theoretical 
interest, which extends to similar problems in other finite many-body systems 
as atomic clusters, quantum dots and atomic gases in traps, the progress in 
this direction would have practical implications for the development of new 
approaches to the solution of the quantum many-body problem. 

Instead of making an attempt to cover recent ideas in a systematic way, 
we show three examples - puzzles which emerge from numerical modeling of 
the nuclear many-body problem in restricted Hilbert space. In all three cases, 
the effect is very clear but its full understanding requires significant theoret¬ 
ical efforts and is not completed until now although below we give plausible 
explanations. 

2 Applying exponential convergence 

The dimensions of shell model spaces increase dramatically with the number 
of single-particle orbitals included. This precludes the full shell model diago- 
nalization in practically interesting regions of the nuclear chart. Even in the 
/p-shell one still awaits for the full calculation for all nuclei. On the other 
hand, such a full solution would provide too much unphysical information 
which is not observable and, moreover, unstable with respect to small varia- 
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Figure 1: Ground state energy of ^®Cr versus the dimension of trucated shell model space 
(total JT-dimension is 232514). 


tions of the interaction hamiltonian that is never known with high accuracy. 
In reality we are interested in the properties of relatively few individual states 
being satisfied with a statistical description for the rest of the spectrum. The 
average properties of excited states were studied in the shell model framework 
with the methods of statistical spectroscopyB and, with direct diagonalization, 
in relation to quantum chaos for atoms 0 and nuclei □. One can expect that 
a reasonable truncation of the shell model space could be possible when the 
influence of the remaining part of the basis is accounted for in some average 
sense. 

Let us consider a behavior of the energy eigenvalue for the ground, or one 
of low-lying, states, as a function of the matrix dimension used for the exact 
diagonalization in the process of the progressive truncation of the original 
huge shell model hamiltonian matrix. The generic picture is shown in Fig. 
1, taken fromQ. After the initial steep decline, the ground state energy of 
^®Cr monotonously descends to the exact value; the convergence is almost 
precisely exponential. This is just one of many available examples, both for 
the realistic shell model and for random Gaussian matrices. The property of 
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exponential convergence seems to be universal □. The full analysis of validity 
can be performed for tridiagonal matrices with a smooth change of matrix 
elements along the diagonal. 

The underlying physics is based on the saturation property of thn^nergy 
dispersion of simple, but spin-isospin {JT) projected, basis states cI’lI. The 
centroid Ek and the width of the basis state |fc) can be found prior to the 
diagonalization in terms of the matrix elements of the shell model hamiltonian: 

E^=Hkk, ( 1 ) 

k^l 

The widths ak are nearly constant for all states of a given JT-class, essentially 
because of the geometric chaoticity of angular momenta couplirig of individual 
particles. This justifies the recipe of statistical spectroscopy □ dealing with 
the centroid E and average width a of each shell model partition. Expanding 
states |fc) in the eigenbasis |a) of the hamiltonian, \k) = 
find the strength function Ek{E) = J2a ~ ^a) which also revealsu 

the saturation property as a function of Ek- The generic shape of the strength 
function evolves, with the interaction strength increasing, from the standard 
Breit-Wigner to the Gaussian one HI . Among various consequences of this 
evolution, one can mention □ the narrowing of the widths of multiple giant 
resonances, r„ —> ^/nTi. In the strong coupling limit, the spreading width 
stabilizes on the level of T 2(t. The strength fragmented to the states at an 
energy di^nce > T should become less and less important. Earlier we have 
suggestedllo the truncation scheme based on this idea. Now we can complement 
this with the exponential extrapolatioiato the exact shell model result. 

The detailed shell model analysis HI has established that the tails of the 
strength function fall off almost pure exponentially as a function of the energy 
distance fropLthe centroid. This phenomenon, which reminds the exponential 
localization til of electronic coordinate wave functions in disordered solids, is 
also related to the nonexponential decay of nonstationary states at the early 
time stage. Based on this result, we can assume that, in the inverted prob¬ 
lem of the composition of the eigenstate |a), the remote basis states \k) give 
exponentially small contributions. This is indeed seen in numerous examples. 

A practical application of the exponential convergence was recently^evel- 
oped for the calculation of ground state energies of nuclei in the /p-shellB. One 
starts with the usual partition of the shell model space and calculates the av¬ 
erage quantities (i- The configurations are ordered according to their energy 
centroids E (this order might be very different from that in the particle-hole 
scheme). After that the diagonalization is performed with progressive inclusion 
of new partitions in their “natural” order. At an energy distance of (3 — 4) S’ 
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from the original centroid, the exponential regime sets inij for low-lying states, 
see Fig. 1, so that extrapolating the energy dependence on the truncated di¬ 
mension n as E(n) = C + B exp(—yn), we can go to the limit oi n = N, the 
full space dimension. 

Using the FPD6 interaction 0, we calculated in this way ground state 
energies, spins and isospins for all the lowest |A(iV — Z)\ nuclides from ^^Ca 
to ®®Ni. Spins and isospins are reproduced correctly, except for the case of 
^®Ti where the three levels with J = 3/2, 5/2 and 7/2 are within 100 keV both 
in the experiment and in our calculation. As usual for such calculations, the 
shell model energies (relative toJdie ^°Ca) require Coulomb corrections E3 and 
additional monopole correctionsEj taking into account the smooth evolution of 
the mean field with the valence particle number. These corrections lead to the 
energy shift without changing the wave functions. The resulting mean square 
deviation of ground state energies from the data is 0.27 MeV. One can conclude 
that the exponential convergence method is a powerful tool to be used in the 
shell model framework for the cases when the full calculation is not feasible. 
The next step in this direction should extend the method to the calculation of 
observables and transition probabilities. 


3 Apparently ordered spectra from random interactions 

The second puzzle was-formulated in the paper 0 that triggered an extensive 
theoretical discussion E3’0. Consider a finite shell-model space with a rota- 
tionally (or/and isospin) invariant two-body interaction. The interaction is 
fully characterized by a set of parameters UL/(jij 2 ; J 3 J 4 ) corresponding to the 
scattering (j 3 ,j 4 ) <-*■ (ji, J 2 ) of the fermion pair in the channel with total spin 
L and isospin I conserved in the process. Let us randomly select these entries 
from a matrix ensemble which is more or less arbitrary but hermitian, real 
and symmetric with respect to the sign of the matrix elements Vli- If, for 
simplicity, the single-particle energies are kept degenerate, what will be the 
distribution function of the quantum numbers of total spin J and isospin T of 
the many-body ground state generated by this ensemble? 

(i) The first idea coming to our mind is that any JT-set has a chance to 
have the lowest energy so that the resulting probability is merely determined 
by the number of available levels with given J and T in the Hilbert space. 
For example, for one kind of particles, according to a traditional consideration 
of the Fermi-gas level density, the total number Af{J) of levels with given J 
can be estimated as A4tat(</) cx (2J -|- 1) exp[—J(J -|- l)/0] where 0 is related 
to the statistical moment inertia determined by the average value of m^, the 
squared single-particle angular momentum projection, in the available space. 
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Figure 2: Fraction fj of ground states with spin J in the uniform ensemble of random two- 
body interactions for = 6 particles on the level j = 15/2; dashed line shows statistical 
multiplicities. 


The maximum of the statistical distribution corresponds to 2 J + 1 = (20)^/^. 

However, this idea turns out to be wrong. As shown in Ref. 0 and con¬ 
firmed by many authors, the ground state spin is predominantly (typically 
with probability exceeding 50%) Jq = 0, although the fraction of states of 
spin J = 0 is usually quite low. The existence of the effect is very robust and 
insensitive to the peculiarities of the ensemble. Its magnitude depends on the 
choice of the ensemble and can exceed 90%. The preponderance of Jq = 0 was 
found also in interacting boson models Ea. In many cases the fraction of the 
ground states with maximum possible spin, Jq = Jm, is also enhanced (the 
statistical fraction of such states is very low; in a single-j model, Fig. 2, the 
state with J = is unique). We know that all even-even nuclei have Jq = 0. 
Is this-Tact originated from pairing forces as suggested by the classics of the 
fieldEjlHj, or will the same pattern appear with nearly any physically allowed 
interaction? 

(ii) Statistical spectroscopy! teaches us to characterize the general features 
of the spectra by the lowest moments of the hamiltonian (centroid, width 
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Figure 3: Widths <jj of the level density pj{E) in the single-j level shell model with random 
two-body interactions. 


and so on). Comparing the statistical widths (t(J) of the subclasses with 
various values of spin J, we may expect that if a class of states with given 
J reveals the largest width (even if the deep reason of that is still unclear) 
the states of this class will be most probably the ground states of the system. 
However, for a majority of ensembles, this conjecture fails. Fig. 3 shows 
the widths aj in the single j-level space for the ensemble of matrix elements 
Vi,L = 0, 2,...,2j — 1, uniformly distributed between —1 and +1 (a system 
of 6 identical particles). Although the ensemble leads to the dominance of 
Jo = 0, the statistical width (Tj=o is not maximum ; in some cases even a 2 > 
CTo- Moreover, to get a signihcant excess of the ground state probability, the 
corresponding width should be considerably greater than others which almost 
never happens. 

If the level density at energy E for spin J is pj{E), the probability of 
finding a state of spin J as a ground state can be formally defined as 



dE 


Xj{E) 


j'^j 


( 2 ) 
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where, in terms of the level densities normalized to 1, gj{E) = pj{E)/JVj, 

Xj{E)=i^J^ dE'gj{E')j . 

For uncorrelated densities pj, eq. (||) indeed prefers the class of states with 
the greatest width. However, this conclusion is not valid since the densities are 
strongly correlated being determined by the same interaction matrix elements. 
The task of calculating the many-point correlation function of level densities 
is very hard. 

(iii) Another alluring idea is that the dominance of Jo = 0 is related to 
the time-reversal invariant character of the random hamiltonian. If so, the 
output could be different for a rotationally invariant, hermitian but complex 
hamiltonian. Physically this can be associated with the fact that the presence 
of Jz ^ 0 acts as if time-reversal symmetry were broken by selecting a sense of 
rotation; the corresponding quasi-Goldstone mode would be rotation restoring 
symmetry by the transformation to another projection J^. This idea in the 
simplest form (introducing an imaginary part of the random matrix element 
Vl) does not workli3 because the ensemble average eliminates all imaginary 
terms along with the odd powers of Vl- Still, the idea is promising if associated 
with the spontaneous symmetry breaking which can be accomplished by the 
consideration of the body-fixed frame, see below. 

(iv) In the first paper td on the subject, see also til, it was suggested that 
usual pairing correlations and the phonon collectivity emerge somehow from 
the random forces. This statement is correct in a limited sense. Indeed, each 
realization of the two-body hamiltonian in a many-level shell model space gen¬ 
erates its own mean field. Then it is possible to construct the superposition 
of the particle-hole operators of a given multipolarity (a generalized phonon) 
which would maximize the coherence and give an enhanced transition prob¬ 
ability. In a simiW way one can look for the specially selected generalized 
seniority operator EH to enhance the pair transfer processes. Those operators 
are different in different copies of the ensemble. A comparison with a standard 
paired state shows that its overlap with the_ground states of random interaii- 
tions is quite small, both in a single j-caseEII and in a realistic shell modeled. 
The phonon collectivity with a fixed multipole operator is also absent. 

(v) Currently the only plausible explanation of the pre fere nce of Jq = 
0 ground states in randomly interacting systems is based Edo on the idea 
of geometric chaoticity. With a random interaction we do not expect any 
specific shape of the mean field to be singled out. The wave functions are very 
complicated combinations of shell model basis states. Therefore a statistical 
approach seems to be suitable which looks for the single-particle density matrix 
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with maximum entropy under constraints of fixed particle number N and total 
spin (isospin). In a single j-model, the density matrix is diagonal for the aligned 
state with the total projection M = J (analog of the body-fixed frame). Its 
eigenvalues give the single-particle occupation numbers = [exp ( 7 m — ^j.) + 
1]“^, where the Lagrange multipliers of chemical potential fi and cranking 
frequency 7 fix the average values of N and M. The expectation value of the 
total hamiltonian calculated with such occupation numbers gives the simple 
approximation for the average yrast line. Depending on the sign of the effective 
moment of inertia for the given set of Vl, this leads to Jq = 0 or Jq = Jm, 
a normal or inverted band, respectively. Thus we come to a trivial geometric 
mechanism of the preference for the edge values of the ground state spin. With 
some improvements, one can reproduce average empirical results. 

The energy values estimated with this statistical approach correlate well 
with ^e exact numerical values although there is a small systematic discrep¬ 
ancy in the probability fj for small J (the agreement is almost perfect for 
high J including Jm)- One source of the deviation is in the approximation of 
the expectation value {fimfim') by the product of two statistical average values 
nmn-m'■ In a given wave function the presence of the mean field makes the 
occupancies slightly dynamically correlated. Apart from that, there is indeed 
some coherence generated apparently by the off-diagonal matrix elements of 
the interaction in higher (even) orders. This brings in a small excess of the 
overlap of the “random” ground state with the paired one cornpared to a pure 
statistical (in the sense of the random matrix theory) estimate cHl. By the same 
reason, the percentage /o of the ground states with Jq — 0 increases when 
going from the single-j case to a more realistic shell model scheme, especially 
to the set of many spin 1/2 levels, when the role of the off-diagonal pair trans¬ 
fer elements is the most important. This can be seen also in the fact of the 
overwhelming percentage of the lowest isospin in the ground states. Another 
example is givenE3 by the specific ensemble which includes only random pairing 
matrix elements; the sign-independent effect of the off-diagonal pair transfers 
leads to the percentage /o > 90%. 

Such dynamical effects are still not fully understood albeit they may be 
the most interesting and essential for many-body physics in finite quantum 
systems. The ideas oLsolving exactly the coherent parts of the interaction 
(for example, pairing cil) and accounting in a statistical spirit for incoherent 
collision-like processes are in the air promising a new interesting development 
in the near future. 



Figure 4: Dynamics of complex energies E—{1/2)T for a system of 3 fermions on 8 equidistant 
orbitals interacting via random interaction; the upper orbital has a single-particle width 7 , 
and the resonances move as 7 increases. 


4 Approaching the continuum 

The third puzzle comes from an attempt to generalize the shell-model approach 
for loosely bound or unbound nuclei where the entire dynamics take place on 
the edge of or already within the continuum. A progress in this direction 
is essential both for nuclear physics far from stability and for astrophysics. 
Since the various versions of the shell model with the discrete spectrum work 
exceedingly well for stable nuclei, it is tempting to consider a realistic-mean 
field where some single-particle orbitals are resonances in the continuumEj and 
include the residual interaction in order to obtain the observable positions and 
widths of the many-body states. 

One example of what happens in such a problem is given by Fig. 4. We 
assume that the single-particle levels have some decay width; for simplicity 
we attribute here a significant width to one upper level shifting its energy to 
the complex plane, e ^ e — ( 1 / 2 ) 7 . Let us switch on a two-body interaction 
with real random matrix elements V (here we do not observe any conservation 
laws so that all pairs of orbitals are mutually coupled) and find the eigenvalues 
of the many-body system. Fig. 4 corresponds to the case of 3 particles on 
8 equidistant orbitals; the trajectories of the complex energies £ = E-(i/2)r 
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are shown as functions of increasing instability 7 of the upper orbital. Instead 
of complete chaos, we see a more and more regular pattern as 7 increases. 
Interestingly enough, 21 energies move almost parallel to each other into the 
complex plane whereas 35 states have a very small width. 

It is easy to understand this dynamics of complex energies (or poles of the 
scattering matrix). The total number of many-body states in this truncated 
space is 8!/(3!4!) = 56. In the limit of large 7 , any state which has the upper 
orbital filled, even with a low probability, will decay very fast. The number of 
such states corresponds to a number of combinations of the remaining particles 
within the rest of space, 7!/(2!5!) = 21. The increasing original decay width is 
distributed over the many-body states, since the imaginary part of the trace 
of the hamiltonian is preserved. As 7 grows, the “self-organization” occurs: 
fast and slow decaying states are separated in time. In a reaction populating 
the system, one would see two distinct time scales, corresponding to direct and 
compound processes. Thus, coupling to the continuum can bring order in a 
system governed by a random hamiltonian. 

The physics we are Icpking at here was extensively discussed earlier from 
a different viewpoint E30’Ea. In that approach one starts with the set of many- 
body states I a) formed by a normal hermitian interaction hamiltonian H. The 
coupling to the continuum is given by the antihermitian part of the effective 
energy-dependent hamiltonian. 



(4) 


C 


Here the sum runs over all decay channels c that are open at a given energy, 
and A% is the decay vertex of an intrinsic state |a) into a channel c. The 



unitarity requirements. The complex eigenvalues of the total effective hamil¬ 
tonian Ti = H — {il2)W give the resonance energies. In the weak continuum 
coupling regime, IT is a perturbation providing narrow resonance widths. As 
this coupling becomes strong, a phase transition occurs to the overlap regime 
with Ericson fluctuations of cross sections and separation of the time scales. 
A number of states (equal to the number of open channels) gives rise to broad 
short-lived resonances absorbing the lion’s share of the total width while the re¬ 
maining compound states become long-lived and reveal internal thermalization 
and equilibration. This phenomenon, being an analog of the Dicke superadi- 
^ce in optics 123, has interesting applications to physics of giant resonances 

In the example above, strictly speaking, all many-body states are nonsta- 
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tionary. In reality there exists a set of threshold energies determined by 
the Q-values of a reaction in channel c. The amplitudes depend on running 
energy E and have a branching point at threshold for example in the case 
of decay into an s-wave in the continuum, (x {E — Therefore, 

in principle, the hybrid approach combining the shell model with the effective 
nonhermitian hamiltonian allows for a self-consistent calculation of discrete 
levels, resonances and reaction cross sections. Being technically very difficult, 
this problem is of vital importance for physics of weakly bound nuclei_Far 
away of thresholds, the method of a nonhermitian hamiltonian was usedEj for 
the microscopic derivation and analysis of the kinetics of resonanc^opulation 
and decay. Two phenomena, the loss of the collective strength S, and the 
restoration of isospin purity at high excitation energy^ naturally follow from 
this consideration. 

5 Conclusion 

Three “puzzles” briefly discussed above show significant gaps in current theory 
of nuclear structure and reactions. The blunt diagonalization of huge shell- 
model matrices cannot be an optimal way of solving the nuclear many-body 
problem. Even randomly taken but geometrically correct interactions generate 
some features of observed regular spectra. The presence of continuum aligns 
the intrinsic states along the new “axis” related to their ability to decay. Those 
are just particular examples of new avenues which should be actively studied. 
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